Search for Dark Matter in the Sun with the ANTARES 
Neutrino Telescope in the CMSSM and mUED 

frameworks 

J. D. Zornoza on behalf of the ANTARES collaboration 

IFIC - Institute* de Fisica Corpuscular, Ed. Instit. de Investigation, AC 22085, 

E-46015, Valencia, Spain 



Abstract 

ANTARES is the first neutrino telescope in the sea. It consists of a three- 
dimensional array of 885 photomultipliers to collect the Cherenkov light in- 
duced by relativistic muons produced in CC interactions of high energy neu- 
trinos. One of the main scientific goals of the experiment is the search for 
dark matter. We present here the analysis of data taken during 2007 and 2008 
to look for a WIMP signal in the Sun. WIMPs are one of the most popular 
scenarios to explain the dark matter content of the Universe. They would 
acumulate in massive objects like the Sun or the Galactic Center and their 
self-annihilation would produce (directly or indirectly) high energy neutrinos 
detectable by neutrino telescopes. Contrary to other indirect searches (like 
with gamma rays or positrons), the search for neutrinos in the Sun is free 
from other astrophysical contributions, so the interpretation of a potential 
signal in terms of dark matter is much more robust. 
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1. Introduction 

One of the most relevant questions in physics is the nature of dark mat- 
ter. Since the first hints about this issue were pointed out by Zwicky in 
1933 [1], a rich set of experimental proofs have been obtained to confirm the 
fact that most of the matter in the Universe is non-luminous. These proofs 
include the observations of WMAP the results on the Big Bang Nucle- 
osynthesis jjjj, the rotation curves of galaxies |I[ and the studies of highly 
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red-shifted la supernovae [5j]. The conclusion is that our Universe is made 
of 73% of dark energy and 27% of matter. Moreover, most of the matter 
(about 80%) is non-baryonic, whose nature is unknown. We know that this 
non-barionic component cannot be dominated by neutrinos, since they travel 
at relativistic velocities and therefore they cannot explain the structure of 
the Universe observed today. Consequently, the candidates for explaining 
most of dark matter are not in the Standard Model. There are several mod- 
els to explain such particles. One of the most accepted candidates come 
from Weakly- Interacting Massive Particles (WIMPs), following the require- 
ments of having an interaction cross-section of the order of that for weak 
interactions (to reproduce the right relic density) and being massive (to ac- 
count of the gravitational effects). They also have to be stable. Among the 
frameworks which are able to provide such particles we will study here two 
of the most popular ones: SuperSymmetry (SUSY) and Universal Extra Di- 
mensions (UED). In particular, we have considered two constrained versions: 
CMSSM p and mUED [7] . The candidate particles in each case are the neu- 
tralino and the lightest Kaluza-Klein particle. Their stability is explained by 
the conservation of the R-parity in the case of SUSY and KK-parity in the 
case of UED. 

The potential detection of WIMPs in neutrino telescopes is based on the 
hypothesis that they would accumulate in massive objects like the Sun, the 
Galactic Center or the Earth due to the loss of energy through elastic scat- 
tering. If they are Majorana particles they would self-annihilate and produce 
secondaries which in turn would yield high energy neutrinos detectable by 
neutrino telescopes. For the case of Kaluza-Klein particles, direct production 
of neutrinos is also possible. In the analysis presented here we focus on the 
Sun. This is very advantageous since a potential detection of high energy 
neutrinos from the Sun direction would be a very clean signal of dark matter. 
In constrast to other indirect searches, there would not be other astrophys- 
ical explantions, as it happens for instance, when an excess of gammas or 
positrons is detected, for which other astrophysical scenarios have to be taken 
into account. 

The analysis presented here is an binned search in the data taken by the 
ANTARES neutrino telescope during 2007 and 2008. 



2 



2. The ANTARES neutrino telescope 



ANTARES [8] is the largest neutrino telescope in the Northern hemi- 
sphere. It is located at (42° 48' N, 6° 10' E) at a depth of 2475 m in the 
Mediterranean Sea, 42 km off-shore Toulon, France. It consists of 885 10-inch 
photomultipliers (PMTs) distributed over 12 vertical lines together forming 
a three-dimensional array. Each PMT is housed in a pressure-resistant glass 
sphere (Optical Module, OM), together with a gel for optical coupling and 
a mu-metal cage to attenuate the Earth's magnetic field. The OMs are ar- 
ranged in triplets (storeys) for a better rejection of the optical background. 
There are 25 storeys per line except for one of the lines which has only 20 
storeys, since the upper part of this line contains devices for acoustic detec- 
tion tests. The PMTs are pointing 45° downward for a better detection of 
upgoing events. The length of the lines is 450 m and the horizontal distance 
between neighbouring lines is 60-75 m. There is also an additional line with 
devices aimed to monitor the environment (sea current velocity, pressure, 
temperature, salinity, etc.) 

The detection principle is as follows. When a high energy neutrino in- 
teracts via a charged current (CC) interaction close or inside the detector 
(in the water or in the rock below), a relativistic muon is produced. This 
muon induces Cherenkov light when travelling through the water, which is 
detected by the photomultipliers. The information of the time [9] and posi- 
tion is used in order to reconstruct the direction of the muon and therefore 
the direction of the neutrino^. The track of the muon is reconstructed with a 
robust fitting algorithm based on a x 2 test which uses the position and time 



of the hits (lOj. It gives an angular resolution of ~2 degrees at energies of 
tens of GeV. The x 2 provided by the fit algorithm is used for the selection 
of well reconstructed tracks, as explained in Section HJ 



3. Data and Monte Carlo simulation 

In this analysis we have used data taken in 2007 and 2008. During 2007, 
there were only 5 lines connected. In May 2008 the detector was completed. 
There were three different configurations: 9, 10 and 12 lines (since May 



1 There are other possible signatures, not used in this analysis: CC interactions of 
electron and tau neutrinos as well as the NC interactions of all flavours, produce shower 
events, which, given the granularitiy of the detector, are seen as point sources of light. 
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2008). The number of active days was 185.5 and 189.8 for 2007 and 2008, 
respectively. The un-filtered data are dominated by the background from 
to atmospheric muons, produced by the interaction of cosmic rays in the 
atmosphere. Most of this background is avoided by accepting only up-going 
events (since muons cannot traverse the Earth) but a fraction of them are 
mis-reconstructed as up going. The other source of background is due to 
atmospheric neutrinos, i.e. neutrinos produced in the interactions of cosmic 
rays in the atmosphere. In this case, the neutrinos can traverse the Earth, 
so it is an irreducible background which is discriminated using the fact that 
this is a diffuse flux, while the signal we are looking for is peaked in the Sun's 
direction. In Fig. [T] we show a comparison between data and simulation of 
these two types of backgrounds. It is important to note that the estimation 
of the background used in the analysis comes directly from scrambled data, 
in order to reduce the effect of systematic uncertainties associated with the 
efficiency of the detector and the assumed flux. 



In order to simulate the signal, the WimpSim package is used [111]. This 



program allows for a computation of the neutrino flux arriving at the Earth 
for different annhiliation channels. The main physical processes are taken 
into account (interactions in the Sun, regeneration of tau leptons, oscillations 
in the propagation, etc.) Once the neutrino flux is computed for the channels 
we are interested in, we can weight the events according to the branching 
ratios of the model under study. As already mentioned, the two frameworks 
we have considered here are CMSSM and mUED. The main channels for 
our detector are Vy + Vy~, bb and rf for CMSSM and cc, bb, rf, ti and vv 
for mUED. In addition to the incoming flux, we have to estimate, through 
simulation, the response of the detector, i.e. the effective area, which relates 
this flux with the number of detected events for a given set of cuts. 



4. Cut optimization 

As mentioned above, we have followed a binned search strategy, i.e. we 
compare the number of events observed in a cone around the Sun direction 
with the expected background, looking for an excess which cannot be ac- 
counted by this background. Since the neutrino annihilation takes place in 
the core of the Sun (which has a size much smaller than the angular resolution 
of the detector), it can be considered as a point source. The optimization 
of the cuts is done following a blinding policy. In our case, we optimize the 
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Figure 1: Comparison between data and simulation: elevation of events. The distributions 
show the simulated atmospheric muons (red), the simulatd atmospheric upgoing- neutrinos 
(blue), and the reconstructed data (black). The ratio of data over simulation is shown 
below the corresponding plot. 
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Figure 2: Normalized effective area for different WIMP masses. 



selection cuts, the search cone size and the quality parameter x 2 i n order to 



minimize the average upper limit 12 



r.90% 

■.in' ;. _ /' /-q 



where fi 90% is the average limit in the number of events, /l eff (M W iMp) is the 
effective area and T e g is the active time. The effective area for the cuts 
Xt < 1-4 and ^ < 3° is shown in Fig. [2j It includes the Sun visibility (i.e. 
how much time the Sun is below the horizon) and the bin efficiency (i.e. how 
much signal is kept in the search cone). 

5. Results 

The average upper limits expected for different models are shown in Fig. [3] 
as a function of the WIMP mass. For the case of CMSSM there is a wide 
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Figure 3: Sensitivity in the neutrino flux as a function of the WIMP mass. (Preliminary). 



spread in the branching ratios for the different models inside the framework. 
Therefore, they are presented separately for the main channels rf , 

bb). The first two channels offer the best limits since the amount of neutrinos 
produced in these channels is higher. The spread for the mUED framework 
is smaller, so we present the sensitivity for one particular realization. In 
this case, the signal is dominated by the rf channel, which explains why we 
obtain similar limits to the corresponding channel in CMSSM. 

These limits can be used to constrain the cross section of interaction 
between neutralinos or KK particles and protons, in particular for the spin 
dependent part. The differential neutrino flux is related to the annihilation 
cross section T by 

dcj) u _ T dN v 

where d is the distance Sun-Earth and dN u /dE u is the differential number 
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of neutrino events for each channel. If we assume that equlibrium between 
capture and annihilation has been reached in the Sun, we can write (for 
definition of C® see Eq. [4J 

r - % (3) 

On the other hand, the capture rate is related to the spin-dependent 
scattering cross-section between a WIMP and a proton as 



C® ~ 3.35 x lO^s" 1 x ( 



Plocal 



0.3 GeV ■cm- 



x 



270km -s-^ ( a H)S D \ x / TeV , . {] 



^local / V 10 6 P b / \ M WIMP 

where pi oca i = 0.3 GeV-cm -3 is the local density of WIMPs, i>i oca i = 270 
km-s" 1 is the local mean velocity of WIMPs assuming a Maxwell-Boltzmann 
velocity distribution, 

M WIMP is the mass of the dark matter particle and (Th,sd is the spin-dependent 
scattering cross-section between a WIMP and a proton. 

We have also used the package SuperBayes [13J and a modified version 



of it [14| in order to scan the available parameter space of the CMSSM and 
mUED models, which are also plotted. It can be seen that neutrino telescopes 
like ANTARES offer the best limits for spin-dependent cross- sect ion. 



6. Conclusions 

We have presented the sensitivity of the search for neutralinos in the Sun 
with the ANTARES telescope. We have used data of 2007 and 2008, when 
the detector configuration had 5 lines (2007) and 10, 9 or 12 (2008). We have 
studied two frameworks, CMSSM and mUED, and the average limits on the 
neutrino flux and the expected limits for the spin dependent cross section of 
the interaction of WIMPs with protons. 
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Figure 4: Experimental limits in the spin-dependent cross-section uh.sd compared to 
the theoretical parameter space allowed by the experimental constraints. The upper plot 
corresponds to the CMSSM model. The bottom plot corresponds to the rnUED model, 
where the light and dark areas correspond to the one and two sigma exclusion regions, 
respectively. (Preliminary). 



